A neutral atom frequency reference in the deep UV with 10 15 range uncertainty 
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We present an assessment of the (6s 2 ) 1 So -O- (6s7s) 3 Po clock transition frequency in 199 Hg with 
an uncertainty reduction of nearly three orders of magnitude and demonstrate an atomic quality 
factor, Q, of ~10 14 . The 199 Hg atoms are confined in a vertical lattice trap with light at the newly 
determined magic wavelength of 362. 5697±0. 0011 nm and at a lattice depth of 20 Er. The atoms are 
loaded from a single stage magneto-optical trap with cooling light at 253.7 nm. The high Q factor 
is obtained with an 80 ms Rabi pulse at 265.6 nm. The frequency of the clock transition is found 
to be 1 128 575 290 808 162.0 ± 6.4 (sys.) ± 0.3 (stat.)Hz (fractional uncertainty = 5.7xl0~ 15 ). 
Neither an atom number nor second order Zeeman dependence have yet to be detected. Only three 
laser wavelengths are used for the cooling, lattice trapping, probing and detection. 

PACS numbers: 32.30.Jc 37.10.Jk 32.10.Dk 42.62.Fi 37.10.Do 



The advance in performance of atomic clocks over re- 
cent decades is an impressive accomplishment. The ac- 
curacy of microwave clocks has improved at a rate of 
~10 per decade and now optical clocks are showing an 
improvement at a rate > 10 2 per decade [l[. Optical 
clock technology is still relatively young; yet four atomic 
species have demonstrated line-center frequency uncer- 
tainties at or below 10 -16 Q-Sj], with the potential of 
several more to follow [BBl • This will enable further pre- 
cision in tests searching for new physics in the low energy 
regime 0-ll|. 

State of the art optical clocks demand an exception- 
ally high line-Q, since the attainable frequency stability 
of a local oscillator locked to an atomic transition is in- 
versely proportional to this quality factor. As Q = v/Av, 
where v is the clock transition frequency and Av is the 
excitation linewidth, it is advantageous to have a car- 
rier frequency as high as practicality allows. In a neutral 
mercury clock this frequency is approximately 1.128 PHz, 
lying in the deep ultraviolet, hence favourable in this re- 
spect. 

Ion clocks have demonstrated impressive results with 
respect to line frequency uncertainties 0, [H, El and 
progress is expected to continue. A prime driver for 
neutral atom clocks is the significantly higher number 
of quantum absorbers used for the atom— probe inter- 
action. This provides an additional measurement lever 
by reducing the integration time required for assessing 
various clock shifts. The tightly bound atoms in an op- 
tical lattice trap become heavily immune to Doppler and 
photon- recoil effects (Lamb-Dicke regime) . In the process 
of constraining the atoms, one needs to shift the upper 
and lower clock state energies by equal amounts by main- 
taining the lattice light at the magic wavelength [ljj . We 
previously reported the measurement of the magic wave- 
length in 199 Hg with 0.21 nm uncertainty @, making it 
the third element to be tested as an optical lattice clock. 
Like other alkaline-earth-metal type elements mercury 
has a doubly forbidden transition between the ground 



(ns 2 ) 1 Sq and the excited (nsnp) 3 Pq levels (n — 6 for 
Hg). It has several favourable characteristics for its use 
as a primary frequency reference; the most significant be- 
ing that its sensitivity to blackbody radiation (BBR) is 
more than an order of magnitude lower than that of Sr 
and Yb [!, HH . The uncertainty due to the BBR dom- 
inates the frequency uncertainty budgets of these two 
clocks [HI, [ljj- Furthermore, the absence of a high tem- 
perature oven for Hg helps to reduce the temperature 
variations in the vicinity of the trapped atoms, thus help- 
ing to reduce the BBR uncertainty further. 

Initial clock transition spectral widths with 199 Hg were 
about 2 kHz due to derealization of atomic states across 
lattice sites 0, fl8| . Here we have increased the lattice 
depth to 20 Er (recoil energy = h 2 kf /2m; ki is the wave 
number; m is the atomic mass) forming more tightly 
bound atoms. This has enabled a linewidth reduction to 
^150 Hz from which a series of light shift measurements 
permitted an improved estimate of the magic wavelength. 
In combination with Zeeman shift measurements this has 
led to further narrowing, reaching 10 Hz. In this letter we 
report results pertaining to both the 150 Hz and ~10 Hz 
wide lines; highlighting the difference where appropri- 
ate. Our narrowest lines correspond to a quality factor 
of Q — 10 14 , which has only previously been obtained 
in a limited number of systems; e.g., 1(3, I3"22|. We 
have conducted measurements related to various system- 
atic shifts, including the first order light shift, Zeeman 
shifts and atom number dependence. 

The 199 Hg atoms are loaded into a ID vertical op- 
tical lattice from a single-stage magneto-optical trap 
(MOT) that uses laser light at 253.7 nm for the cooling 
I 1 So o 3 Pi). The MOT and lattice trap are described 
in detail in }18| . Some modifications have been made to 
the lattice cavity and to the resonant cavity generating 
the 253.7 nm ultraviolet light, which we outline here. A 
scaled drawing of the apparatus is shown in Fig.JIJa). On 
the left hand side appears the combined MOT chamber 
and optical lattice cavity. The ports in the horizontal 
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and 45° directions provide access for the MOT beams 
(from beneath in the case of the 45° ports). On the right 
hand side is the mercury source chamber, which includes 
a 2D MOT set-up. The 2D MOT is not employed in the 
measurements described here, but its use is anticipated 
in future. Several grams of mercury are maintained at 
about — 40 °C using a dual Peltier stage. Between 50 and 
70 mW of 253.7 nm light is generated from a frequency 
quadrupling scheme for use in the 3D MOT (with the 
variation in power occurring over weeks or months). 
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FIG. 1: (Color online) (a) Drawing of the vacuum system showing 
the 3D MOT (left side) and 2D MOT chambers. Also visible is 
the upper mirror of the optical lattice cavity above the 3D MOT 
chamber; fl., fiuoresence. (b) Timing sequence for the cooling light, 
probe light and 3D MOT magnetic field gradient. The lattice light 
remains on continuously throughout the cycle. Af ; = 80 ms and 
T p = 50 ms for most measurements here; Atj = 9 ms. 



The vertically orientated lattice cavity is comprised of 
two spherical mirrors, both with a radius of curvature of 
250 mm and has a finesse of 210 at 362 nm. The waist 
size (rad.) is 120 /im where the lattice light overlaps the 
MOT cloud and produces a maximum lattice depth of 
25 En (or 9.2 /iK). Light at or near the magic wavelength 
is coupled into the cavity from below and the transmit- 
ted light is used to form a side-lock to maintain constant 
intensity. A tunable Ti:sapphire laser has its output fre- 
quency doubled in a LIB3O5 crystal based resonant cav- 
ity to produce the lattice light. In our previous lattice 
cavity design there was a 45° reflector that formed an re- 
shaped cavity to lift polarization degeneracy fl8| . This 



optic was found to degrade rather quickly under vacuum 
with the incidence of 362 nm light at high power and has 
been removed. 

The 1 S$ <H> 3 Po clock transition lies at 265.6 nm, so 
as in the case for the 253.7 nm light generation, two fre- 
quency doubling stages are employed (see Fig. (2fa) for 
the relevant electronic transitions). The infrared-light 
(IR) is sourced from a distributed feedback semiconduc- 
tor laser, injection locked with light from a fiber laser 
tightly locked to an ultrastable optical cavity [23| . About 
1 mW of 265.6 nm light is produced by the frequency 
quadrupling scheme. For rapid control of the 265.6 nm 
probe light level we employ an AOM, which positively 
shifts the frequency of the light by 180 MHz. Swee pin g 
the frequency of the 265.6 nm light is described in [18]. 
An AOM that is used to suppress noise in the fiber link 
between the ultrastable laser and the main Hg appara- 
tus is used to tune the clock probe frequency. Despite 
the drift rate of the ultrastable laser remaining below 
+30mHzs~ 1 over the last five months, we still find it 
helpful to include a de-drift scheme to keep track of the 
clock transition. This is performed using a direct digital 
synthesizer (DDS) that steers the frequency of the AOM 
in the 1062 nm path. 

The profile of the clock transition is made via detec- 
tion of atoms in the ground state only and with the 
timing sequence shown in Fig. [IJb). A broad scan of 
the transition spectrum is shown in Fig. [21b), where 
the magnetic bias field is made small enough that the 
^Sq) m F = ±1/2 o ( 3 P ) m F = ±1/2 Zeeman com- 
ponents overlap one another. The frequency is offset by 
the value reported in [24| for the 199 Hg transition fre- 
quency; i.e., v c = 1128575290808400Hz. We applied 
1.5 /iW of 265. 6nm light with a e -2 beam radius of 
310 /im (intensity = lOWnr 2 ). The full-width half- 
maximum (FWHM) is 140 Hz and the contrast ^32%. 
We show below that a much narrower transition lies at 
the center of this spectrum. We will henceforth refer to 
this broad profile as the pedestal. This pedestal may 
be an indication that the transverse confinement of the 
atoms should be improved. When applying a dc B-field, 
this line profile separates into two Zeeman components 
[Fig. [He)] with slightly lower contrast and FWHM equal 
to 120 Hz. The Zeeman line separation versus bias field 
strength is displayed in Fig. &[d) (circles) exhibiting a 
slope of 11 Hz /iT -1 . Since the 265. 6nm probe light co- 
propagates with the lattice light, its polarization lies or- 
thogonal to the axial direction of the lattice (i.e. in the 
horizontal plane). The bias _B-field is applied in the in- 
direction seen in Figflja). The probe light polarization 
was confirmed to be mostly linear, but shares a compo- 
nent in both x and y-directions, allowing the possibility 
for both 7r and a transitions. The HHz/^T -1 depen- 
dence corresponds closely to that expected for a transi- 
tions [25| . Ideally the probe light polarization should be 
in the ^-direction, but in the present experiment the ul- 
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tranarrow clock line strength (below) is optimised when 
the polarization has shared x and y components. This is 
still a matter for investigation. 



(6s6p) 



(6s7 s) 3 Si 




B = 0mT 

0.64, , , , , , 

800 -400 400 
Frequency - v c (Hz) 



(a) 



(b) 

0.0 0.2 



c 1.00- 
o , 

| 0.95- 
8 0.90- 

OS 

■S 0.85- 

c 

0.80- 



0.75 1 



-800 



B = 35 uT 

-400 400 
Frequency -v c (Hz) 

(C) 



600 

N 

£ 400 

c 

I 200 
& 

if) 

1-200 

<D 
<D 

M -400 











■ 






Hi 








, M 












- 


i 




♦ It 

• pedestal 



0.4 (G) 
60 

40 

20 

- 

-20 

40 



-20 







20 



40 



Magnetic field strength (^T) 

(d) 



FIG. 2: (Color online) (a) Partial level scheme for 199 Hg with hy- 
perfine splitting of the ground and excited states. The 265.6 nm 
radiation is used to probe the 1 So <-)> 3 i 3 o clock transition, while 
253.7 nm radiation is used for cooling and detection, (b) A spec- 
trum of the 199 Hg clock transition "pedestal" with mp = ±1/2 — > 
mf = ±1/2 Zeeman components overlapped. Ar J =362.573 nm and 
FWHM = 140 Hz. (c) Line profile of the 199 Hg pedestal show- 
ing the separated Zeeman components with FWHM =120 Hz. (d) 
Frequency separation of the Zeeman components for the ultranar- 
row 7r transition and pedestal versus bias B-field. The slope is 
3.1 ±0.3 Hz //T- 1 and 11.1 ± l.THzjiT -1 , respectively. 



Information about the confinement of the atoms in lat- 
tice trap is garnered by examining the sideband spec- 
tra of the Lamb-Dicke spectrum. Although not shown 
here, when we curve- fit to the blue sideband using the 
approach outlined in 26, 27] with the only free parame- 
ters being the lattice depth and the temperature in the 
transverse direction, we find U = 18 Er (or 6.5 /j,K) and 
T r = 7/iK, with associated axial and transverse frequen- 
cies of 64 kHz and 43 Hz, respectively. We estimate the 
temperature of the atoms in the axial direction based on 
the ratio of the blue and red sidebands to be ~4 /iK (with 
P n /P n+ i ~ 2.4). 

With the lattice light set at the magic wavelength (dis- 
cussion below) and a higher resolution scan made across 
the center of the pedestal, we find a much narrower spec- 
tral line. Figures [3ja) and (b) show spectra obtained 
with T p = 50 ms and 80 ms probe pulse durations taken 



with approximately 1 ^W and 0.5 ^W of 265.6 nm probe 
light, respectively. The solid lines are the modelled Rabi 
spectra with il.T p ~ 1.47rrad for both. In Fig.^a) £l.T p 



was chosen in an attempt to match to the sidebands ob- 
served. While the width of the theoretical trace matches 
the data, the frequency of the sidebands in general does 
not, as seen in[3][b). From a measurement of the probe in- 
tensity we estimate the Rabi angle to be il.T p ~ 4.57rrad. 
Despite this, we suspect that the sidebands are of a tech- 
nical origin and are not due to an overdriven Rabi pulse, 
since reducing the probe power does not decrease the 
size of the sidebands. For most line-center measurements 
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FIG. 3: (Color online) Ground state fraction versus probe laser 
detuning frequency for: (a) 50 ms (Q.T P = 1.457rrad) and (b) 80 ms 
(CI. T p = 1.367rrad) Rabi pulses. B=0 and A^RiA m for (a) and (b). 



here we use a 50 ms probe. With an applied bias _B-field 
the narrow transition separates into two Zeeman com- 
ponents. The 1 st order Zeeman dependence is seen in 
Fig.^d) with a slope of 3.1 ± 0.3 Hz fiT^ 1 . The ratio of 
the -B-field dependencies for the pedestal versus the nar- 
row 7r transition is ~3.6, which is close to the ct/tt ratio 
of 3.3 that one expects when calculated from the differ- 
ence in the Lande g-factors of the ground and excited 
states [25| . 

For assessments of the ac Stark shift, the center fre- 
quency of the clock transition is measured at a series 
of lattice wavelengths and lattice depths. The results 
are summarized in Fig. [4ja) . The main graph shows 
data taken with the 10-15 Hz wide spectral lines, while 
the inset shows lower resolution data obtained with the 
pedestal. From these data we determine the magic fre- 
quency to be 826.8546 ± 0.0024 THz (A m = 362.5697 ± 
0.0011 nm). The uncertainty also incorporates the accu- 
racy of the wavemeter used for the measurements. From 
the line fit we find the strength of the light shift is 
(-5.1 ± 0.9) xlO -17 ^ 1 GHz -1 . Another representation 
of the light shift is shown in Fig. SJb) . Here the line cen- 
ter frequency is plotted as a function of the lattice wave- 
length while some variation was applied to the lattice 
depth as indicated by the color scale. The small vari- 
ation in center frequency near 362.570 nm is consistent 
with the above magic wavelength determination. Here 
the strength of the light shift is exlO" 17 ^ 1 GHz" 1 . 

A preliminary assessment of other systematic shifts af- 
fecting the 199 Hg clock transition has been conducted 
and is summarized in Table HI For a temperature uncer- 
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FIG. 4: (Color online) (a) Differential light shift versus lattice 
frequency for the 199 Hg clock transition. The inset shows mea- 
surements made over a larger frequency range. The Stark-shift 
free frequency (wavelength) is 826.8546 ± 0.0024 THz (362.5697 ± 
0.0011 nm) and the slope is -57 mHz .E^ 1 GHz -1 , (b) Clock tran- 
sition frequency versus lattice wavelength with changes of lattice 
depth. 
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tainty of 2 K of the chamber surrounding the atoms, and 
a 100 % uncertainty for the calculated BBR coefficient 
the shift at 290 K is -0.17 ± 0.20 Hz. A measurement of 
line center frequency versus Zeeman component separa- 
tion showed no variation within a statistical uncertainty 
of 1.6 Hz. This is expected since the predicted 2 nd order 
Zeeman dependence is 2.16xl0~ 9 T~ 2 (e.g., B = 10//T 
produces a shift of 0.24 mHz) 

The frequency dependence on atom number has been 
tested by two means: (i) by varying the loading time of 
the MOT and (ii) by varying the level of 253.7 nm light 
used for cooling. Both methods change the number of 
atoms loaded into the lattice trap (and also the density 
since the size of the MOT cloud changes by less than 
5% for the range considered here). The nominal atom 
number is Nq = 2.5 xlO 3 (close to the maximum that 
we achieve). Fig. [5ja) shows the line-center frequency, 
^Hgj versus the relative atom number (N/Nq) obtained by 
varying the MOT loading time from 0.8 to 2.5 s. Fig.JSfb) 
shows z/Hg for a range of cooling light intensities with a 
slope equal to —0.11 ± 0.08 Hz per unit so (s = I/Iq, 
where Jo=102W m -2 ). The two results suggest a non- 
significant density shift at the present resolution. For 
daily assessments of line-center frequency (described be- 
low) there are between 10 and 30 spectra recorded. From 
a plot of center frequency versus relative atom number 
we determine the frequency shift at Nq, which we show 
in Fig. [5fc) for each measurement day where there was 
sufficient variation in N/Nq. From the weighted mean 
the density shift is constrained to 0.26 ± 1.9 Hz, consis- 
tent with the previous two results. Note, the number 
of atoms per site is of order one and in the MOT peak 
atomic density is ~ 6xl0 10 cm~ 3 . 

An ac Stark shift may also arise due to the probe light. 
In Fig. E^d) we show the clock transition frequency ver- 
sus probe power. At our nominal power level of 1 /iW 
the shift lies within the uncertainty of the measurement. 



FIG. 5: Line center frequency versus: (a) relative atom number, 
where 7Vo = 2.5xl0 3 and (b) MOT cooling light intensity, so. (c) 
Frequency shift at jVrj versus measurement day. (d) Line center 
frequency versus probe power. 



TABLE I: Corrections and uncertainties for the 1 Hg clock 
transition 



Effect 


Correction (Hz] 


Uncertainty (Hz) 


Blackbody radiation 


-0.17 


0.2 


1st order Zeeman 


< 1 


0.5 


2nd order Zeeman 


< 1 


1.6 


Scalar light shift (lattice) < 1 


5.9 


Probe light 


0.02 


0.3 


Atom number density 


0.4 


1.9 


Total 


0.2 


6.4 



We also calculate that the probe light shift is in the 
1 — 10 mHz range. Tensor Stark and hyperpolarizability 
shifts due to the lattice light are omitted from the table 
as they are expected to be at least an order of magnitude 
.smaller. 

Frequency comb accuracy has been verified by com- 
paring measurements from two frequency combs (a 
Ti:sapphire comb and a Er:fiber based comb), which 
are steered by a low noise H-maser reference. The two 
combs show agreement to below 0.01 Hz. The H-maser 
reference is continually compared to the LNE-SYRTE 
frequency standard, which has a relative uncertainty of 
2.7x 10" 16 [HI. An example of the reproducibility of the 
clock transition frequency is shown in Fig. [5J where mea- 
surements have been made over a period of three months. 
Each point is the mean line center frequency produced 
from between 10 and 30 spectra recorded on each day 
represented. The measurements were mostly performed 
with the two (7r-transition) Zeeman components over- 
lapped due to the low S/N in the present experiment. 
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FIG. 6: (Color online) Clock transition frequency measured over a 
three month period with respect to the SYRTE primary frequency 
standard for A m = 362.5697 nm. The weighted mean is v c — 238. 2± 
0.3Hz. v c = 1128575290808400Hz. 



The variations in frequency reduce after MJD = 55825 
when we began controlling the lattice light frequency to 
within 300 MHz. The uncertainty of the weighted mean 
of this series is 0.3Hz (a/vn g = 2.5 x 10~ 16 ). Accounting 
for the systematic shifts (Table HJ we evaluate the 199 Hg 
1 S <-> 3 P transition frequency to be 1 128 575 290 808 
162.0 ± 6.4 (sys.) ± 0.3 (stat.) Hz, where the combined 
fractional uncertainty is 5.7xl0 -15 . 

This work establishes the potential of 199 Hg as a high- 
accuracy clock. There remain various means for fur- 
ther gains in accuracy; e.g., the S/N of the clock sig- 
nal should improve as techniques for sideband cooling, 
transverse cooling and atom number normalization are 
implemented. One also expects that a reliable means of 
producing > 150 mW of 253.7 nm radiation will be found 
in which 2D-MOT loading should increase atom numbers 

10-fold nu. 
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